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We theoretically investigate resonant tunneling through S- and U-shaped nanostructured graphene 
nanoribbons. A rich structure of resonant tunneling peaks are found eminating from different quasi- 
bound states in the middle region. The tunneling current can be turned on and off by varying the 
Fermi energy. Tunability of resonant tunneling is realized by changing the width of the left and/or 
right leads and without the use of any external gates. 



PACS numbers: 73.23.-b, 73.40.Gk, 73.40.Sx, 85.30.De 



INTRODUCTION 



Graphene is a single layer of carbon atoms arranged in 
a hexagonal lattice that is the building block for graphite 
material. Recently, graphene samples have been fab- 
ricated experimentally by micro-mechanical cleavage of 
graphite 1 . This material has aroused increasing atten- 
tion due to its novel transport properties that arises 
from its unique band structure: the conduction and 
valence bands meet conically at the two nonequivalent 
Dirac points, called K and K valleys of the Brillouin 
zone, which show opposite chirality. Around the two 
Dirac points, the energy dispersion is linear in momen- 
tum space and is well described by the massless Dirac- 
Weyl equation. The unique energy dispersion leads to 
a number of unusual electronic transport properties in 
graphene such as integer quantum Hall effect at room 
temperature^, finite minimal conductivity, and special 
Andreev reflection^. In graphene electrons can pass 
through potential barriers without any reflection, i.e., 
the Klein tunneling, in contrast to conventional electrons 
that exhibit an exponential decreasing transmission with 
increasing barrier height and/or width. Klein tunnel- 
ing makes it impossible to confine the carriers in single 
layer graphene using the electric gates that are often used 
in the conventional semiconductor two-dimensional elec- 
tron gas (2DEG). Tunneling through single barriers^ and 
double barriers 6 - in graphene has been investigated theo- 
retically. Recent theoretical work demonstrated that the 
magnetic barriers can be used to confine the massless 
Dirac fermion 7 and block the Klein tunneling and dis- 
play a huge magnetoresistance^. Klein tunneling leads 
to a poor rectification effect and the absence of resonant 
tunneling for normal incidence in graphene and limits the 
performance of graphene-based electronic devices. 

In graphene nanoribbons (GNRs), the presence of 
edges can change the energy spectrum of the 7r-electrons 
dramatically. GNRs can be fabricated using conventional 
lithography and etching technique a 10 ^ 1 . The electronic 
properties of a GNR depend very sensitively on the sizes 
and shape of the edges, i.e., zigzag- and armchair-edged 



GNR. The zigzag-edged graphene nanoribbons (ZGNRs) 
and armchair-edged graphene nanoribbons (AGNRs) ex- 
hibit very different band structures. The ZGNRs al- 
ways have localized states appearing at the edge near 
the Dirac point, and therefore exhibit metallic-like be- 
havior. The AGNRs show alternating metallic-like and 
semiconductor-like features alternatively as the width of 
nanoribbons increases^. Such features are very different 
from the conventional semiconductor quantum wire and 
provide a unique way to tailor the transport properties of 
GNRs±^^*±&±L±£. Recently it is found that single elec- 
tronic potential barrier can block the electron tunneling 
in the ZGNRs and confine electrons in between double 
barriers^. The magnetic transport and spatial distribu- 
tion of currents was investigated for unipolar and bipo- 
lar ZGNRs 16 . These previous works have demonstrate 
that the resonant tunneling can be realized using elec- 
tric and/or magnetic barrier. Here, we will address the 
question whether it is possible to realize resonant tunnel- 
ing using only the geometry of the graphene nanoribbons. 
We focus on the S- and U-shaped structures consisting of 
zigzag leads and armchair nanoribbons in between them, 
because the zigzag nanoribbon is always metallic, while 
the bandgap of an armchair nanoribbon depends sensi- 
tively on the width of nanoribbon. This combination of 
zigzag and armchair nanoribbon could lead to interesting 
transport property. 

In this work, we theoretically investigate the transport 
property through U-shaped nanoribbon structures sand- 
wiched by two semi-infinite ZGNRs, i.e., the geometry of 
the graphene nanoribbon in the absence of any external 
nanostructured electric gates.. We will demonstrate res- 
onant tunneling behavior through these graphene nanos- 
tructures whenever the AGNR in the middle part is 
metal-like or semiconductor-like. In those structures, the 
middle structure acts as a barrier between ZGNR leads, 
and an interesting resonant tunneling behavior can be 
seen with varying lengths of the middle part of the U- 
shaped ZGNR or AGNR. 

This paper is organized as follows. In Sec. II, the the- 
oretical model and calculation method are given. In Sec. 
Ill, we present the numerical results and our discussions. 
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FIG. 1: (Color online ) Schematic diagram of the multi-moded 
transmission. 



Finally we give our conclusion in Sec. IV. 



II. MODEL AND CALCULATION METHOD 

To describe the electron transport through the S- or 
U-like structure sandwiched by two semi- infinite ZGNR 
leads, the tight-binding Hamiltonian is adopted 

H = ~ l E fa + H - c ) - l ' E ( c ta + H - c ) > 



<hj> 



(i) 

where ct^ (c^j)) creates (annihilates) an electron on site 
Ri(j), t is the nearest neighbor (< z, j >) hopping en- 
ergy (t & 3.03 eV), and t' is the next-nearest neigh- 
bour (<< z, j >>) hopping energy 19 . The conductance 
of the system is evaluated using the Landauer-Biittiker 
formula 20 , 

G (^) = 7 E t m (^) - ^ (^) = E i w (^) i 2 ' 

(2) 

where is the transmission coefficient from the 

^-th channel in the left lead to the ja-th channel in the 
right lead at the Fermi energy Ep (as shown in Fig. [1]), 
calculated by a recursive Green's function method 21 , 
which is given as 



G s = 



1 



1 



E-H E-Hc-^i 



(3) 



where E is the Fermi energy, H is the effective Hamil- 
tonian that equals the sum of the Hamiltonian of the 
central rectangular ring, He, and the self energies of the 
semi-infinite leads in the left and right leads, E/, + E#. 
The key point of the recursion Green's function method 
is that the middle conductor region is divided into many 
meshes, and consequently reduce the dimensionality of 
the Green's function greatly. Therefore the computation 
of the conductance becomes more pratical. 

The Fermi energy Ep can be tuned experimentally 
through electric top- and/or back-gates 22 . In our cal- 
culations, all physical quantities are dimensionless, e.g., 
the energy Ep and the conductances are in units of t and 



e 2 / (it), respectively. The next-nearest neighbour hop- 
ping term is ignored because it is much weaker than that 
of the nearest neighbor hopping term, i.e., t' = t/10 4 . 

For completeness, here we repeat the essential steps 
in T. Ando's work 21 to obtain recursive Green's function 
Gs • The equation of motion for an ideal lead can be 
written as: 



(E - H )Cj + TC,_i + T T C J+1 = 0, 



(4) 



where Ho is the hamiltonian of a mesh in the ideal lead 
containing M lattice sites, Cj is a vector describing the 
amplitudes of the jth mesh, T and T T are the matrix 
blocks adjacent to H$ in the total tight-binding Hamil- 
tonian matrix. To obtain linearly solutions for Eq. 0J we 
first set Cj = VCq. Substituting this in to Eq. [4] we 
have: 



AC i = (T T ) (Ho — E) Cj — (T T ) TCj.i 
leads to the following eigenvalue problem: 

(T T y X (H -E) -(T T y l T 







(5) 



Cj-i 
(6) 



We get the eigenvalues and corresponding eigenvectors, 
including M right-going and M left-going waves. We de- 
fine the matrix for the eigenvalues and eigenvectors: 



I7(±) = («i(±),-",u M (±)); 



A(±) 



/ Ai (±) 

V 



(7) 



(8) 



A M (±) 

Now, we get the relations between j—th and j —th mesh: 
C i (±) = F(±)^ / C i /(±), (9) 

with F{±) = U (±) A (±) U~ 1 (±). 

Next we consider the scattering in the middle conduc- 
tor (from mesh lto mash N) to both sides of witch an 
ideal wire is attached. By using Eq. 0] and Eq. [9l we 
obtain the effective hamiltonians (note that they are not 
hermitian) of the last mesh in left lead and first mesh in 
right lead, both adjacent to the middle conductor: 

H = H -TF- 1 (-), (10) 

H N+1 =H -T T F(+). (11) 

We can get the total effective hamiltonian matrix: 



f Ho 
-T 



H 



\ 








_ T T 

Hi 









Hn 
-T 






_ T T 

Hn+i 



\ 



(12) 



where Hj = Hj for j = 1, 2, iV in the central region. 
Proceed by using Eq . [3] directly, we can obtain the de- 
sired recursive Green's function Gs • 
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III. RESULTS AND DISCUSSIONS 

First, we discuss the transport property through the 
double-bended (S-shaped) GNR structure as shown in 
Fig. [2] (a) for reference purposes. Fig. [2] (c) shows the 
conductance as a function of the Fermi energy for differ- 
ent lengths of the AGNRs. There are the strong peaks 
in the conductance spectrum corresponding to the reso- 
nant tunneling process through the double-bended GNR 
structure. The underlying physics of such resonant tun- 
neling process is that there are quasi-bound states in the 
middle conductor region which consists of the armchair 
and zigzag nanoribbons, the resonant tunneling occurs 
when the Fermi energy is equal to the energies of the 
quasi-bound states. The quasi-bound states arise from 
the propagation of electrons back and forth in the arm- 
chair region between the zigzag leads. The conductance 
of the perfect ZGNR displays a step-like feature, which 
corresponds to the opening of the different channels as 
the Fermi energy increases (see the green line in Fig. 
Efc)). Note that we did not introduce any external elec- 
tric field and/or gate to generate the barrier in the middle 
part. Therefore, this feature arises from the geometry ef- 
fect, i.e., the different edges of the ZGNR leads and the 
AGNR conductor. As the Fermi energy increases, more 
incident channels are opened, the effect of the edge state 
becomes weaker, and the ratio between peaks and valleys 
of the conductance decreases (see the black and red lines 
in Fig. (21(c)). When the Fermi energy increases further 
and even becomes higher than the bottom of the second 
subband of the middle AGNR conductor (see the black 
vertical line in Fig. [2] (b)), the oscillation of the conduc- 
tance becomes more significant. Fig. [2](c) indicates that 
there are more tunneling peaks as the AGNR length in- 
creases. The incident electron can be completely reflected 
when the incident energy is less than the bottom of the 
second subband of the left ZGNR lead corresponding to 
the case of single mode incidence. 

To obtain a clear physical picture, we plot the distri- 
bution of electron states at specific Fermi energies Ep 
corresponding to the points A and B marked in Fig. [2] 
(c). Figure [2] shows the density distribution of the quasi- 
bound state in the middle AGNR part that contributes 
to the resonant tunneling (see the point A in Fig. [2] (c)). 
The electrons mainly localize at the edges of the mid- 
dle AGNR, at the corner between the AGNR and ZGNR 
leads, and on the edge of AGNR. This quasi-bound state 
is quite different from the electron states of a perfect 
AGNR where no edge state appears due to the two dif- 
ferent kinds of carbon atoms at the edges. But for the 
quasi-bound state, the electron state strongly localizes at 
the edges and especially at one kind of carbon atoms (see 
Figs. [2] (d) and (e)). This special quasi-bound state in 
the middle part contributes dominantly to the resonant 
tunneling process. We also plot the density distribution 
of the quasi-bound state corresponding to the completely 
reflection case (see point B in Fig. Efc)) in Fig. [2](e). For 
this case, the electron state localizes at the edge of the left 



semi-infinite ZGNR lead near the interface between the 
lead and the middle conductor. In this case, the AGNR 
conductor can be used as an electron reflector. For the 
usual graphene nanoribbon, there is no quasi-bound state 
in the middle region since no armchair nanoribbon lo- 
cates in between the zigzag leads. Therefore we can only 
observe the step-like feature of the conductance which 
corresponds to the opening of the channels with higher 
energies, (see the green line in Fig. 3(b)) 

Next, we show how the width of the right lead affects 
the resonant tunneling behavior in Fig. [3l We fix the 
width of the left lead at Np = 20 while tuning the width 
of the right lead. In Fig. O the black solid, the black 
dashed, and the red dashed lines correspond to Nr = 20, 
10, and 40, respectively. Figure [3] (a) demonstates that 
the magnitude of the tunneling peaks decreases whether 
the width of the right lead is wider or narrower than that 
of the left lead. Thus the electrons cannot be fully trans- 
mitted through the nanostructure. The perfect transmis- 
sion only takes place when the right lead has the same 
width as the left lead. When the width of the right lead is 
narrower than that of the left lead (see the black-dashed 
line in Fig. [3]), the first tunneling peak is weakened heav- 
ily. When the Fermi energy increases, the magnitude and 
width of the tunneling peaks increase gradually. These 
increases are caused by the enhanced coupling between 
the quasi-bound states in the middle part and the leads. 
Comparing the case with equal widths of the left and 
right leads shown in Fig. [3] (a), all tunneling peaks for 
the wider (narrower) right leads shift to higher (lower) 
energy, and become sparse as the Fermi energy increases. 
For the wider right lead case, the conductance shows dif- 
ferent features. For the Fermi energy below the bottom 
of the second subband in the right lead (see the first ver- 
tical green line), the variation of the magnitude of the 
tunneling peaks is similar to the case with the narrower 
right lead. When the Fermi energy is above the bottom 
of the second subband in the right lead, the magnitude 
of the tunneling peaks increases abruptly, even approach- 
ing to 1, and becomes wider (see the red dashed peaks 
in Fig. 03(b)). Comparing the conductances through the 
middle AGNR with different lengths (N H ) in Fig. C3 
we find that the tunneling peaks become dense as the 
length N H increases. For longer AGNR (N H = 160), the 
variation of the width of the right lead induces changes 
in the conductance similar to that of the shorter AGNR 
(iV^ = 80). 

Now we turn to consider the U-shaped structure com- 
posed of two double-bend structures as shown in Fig. [H 
Figure [5] (a) shows the contour plot of the conductance 
as a function of the length of the middle ZGNR bridge 
(Ntc) an d the Fermi energy (Ep) for the U-shaped 
structure. From this figure, the tunneling peaks show a 
periodic behavior as a function of the length Ntc- This 
period is exactly the same as that between the met al- 
and semiconductor-like AGNRs. The resonant tunneling 
processes can be divided into two kinds of processes ac- 
cording to the behavior of tunneling peaks varying with 
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FIG. 2: (Color online) (a) Schematic diagram of the double-bended graphene nanoribbon structure, (b) The band structure 
of the perfect armchair graphene nanoribbon (AGNR) with the same middle conductor width (N c ). (c) The conductance as 
a function of the Fermi energy through the double-bended structure, the black-solid and red-dashed lines correspond to the 
different conductor lengths Nh=80 and 160, respectively, the green line is for the perfect zigzag graphene nanoribbon. (d) and 
(e) show the density distributions of quasi-bound electron states marked by the points A and B in Fig. [2] (c) for Nh = 80 
double-bended structure. Nl=Nr=20 and Nc=S. 
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FIG. 3: (Color online) (a) and (b) The conductance as a func- 
tion of the Fermi energy for Nh = 80 and 160, respectively. 
The different lines correspond to the different widths of the 
right lead for fixed width of the left lead A^=20. The blue 
vertical lines correspond to the bottom of the second and third 
subbands for the perfect ZGNR with the same width of the 
right lead A^=40. 



the length Ntc- The first kind of tunneling peaks is 
strong and does not shift when the length of the middle 
bridge increases, while the second kind shows the oppo- 
site behavior, i.e., shifting to lower and/or higher ener- 
gies. Comparing the conductance of the U-shaped struc- 
ture with that of one double-bended (S-shaped) struc- 
ture, the four peaks in Fig. 03(b) correspond exactly to 
the tunneling peaks of the double-bended structure (S- 




Left lead 



Right lead 



FIG. 4: (Color online) Schematic diagram of the U-shaped 
graphene nanoribbon structures which is composed of two 
double-bended nanoribbons. The width of the left (right) 
ZGNR lead is denoted by N L ( R ^ in the unit of the number 
of hexagonal units. The width and length of the left (right) 
AGNR are denoted by N LC (N RC ) and N LH (N RH ). N M 
and Ntc correspond to the width and length of the bridge 
GNR (the black frame) connecting the left AGNR and the 
right AGNR. 
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FIG. 5: (Color online) (a) The contour plot of the conduc- 
tance as function of the length of the ZGNR (Ntc) and 
Fermi energy (Ef) for the U-shaped structure. N L ( R ^ = 20, 
Nlc N rc = 8, and N LH N RH 80. (b) the con- 
ductance as a function of the Fermi energy for the double- 
blended (S-shaped) GNR structure. N L(R) = 20, N c = 8, 
and Nh = 80. The red line denotes the conductance as a 
function of the Fermi energy for the U-shaped structure at 
Ntc = 42. 



shaped) (see Figs. [5] (a) and (b)). Therefore these reso- 
nant tunneling processes are caused by the quasi-bound 
states localized in the two AGNR regions in the middle 
conductor parts, which behave like the resonant tunnel- 
ing through a single barrier. This can be understood 
from the spatial distribution of the quasi-bound states 
contributing to the resonant tunneling process marked by 
point A in Fig. [5](a). The quasi-bound electron state lo- 
calizes at the left and right AGNR of the middle conduc- 
tor. From the inset of Fig. 03(b), one can see that the first 
kind resonant peak near E=0.05t split into two peaks due 
to the coupling between the two lowest localized states in 
the left and right AGNRs through the ZGNR bridge, i.e, 
the bonding and antibonding splitting. Therefore the res- 
onant peaks split into two peaks and show an anticross- 
ing behavior with increasing the length of ZGNR bridge 
that describes the coupling strength. The other resonant 




FIG. 6: (Color online) (a), (b), and (c) show the density 
distributions of the quasi-bound electron states corresponding 
to the points A, B, and C in Fig. [5] for the N c = 44 U-type 
structure. 



tunneling peaks are shaper and are located between the 
strong tunneling peaks mentioned above and vary as the 
length Ntc increases (see Fig. [5] (a)). Those peaks is 
the result of the resonant tunneling processes through 
the double barriers, since the corresponding quasi-bound 
states localize in the ZGNR bridge of the middle conduc- 
tor. These quasi-bound states mainly distribute in the 
ZGNR bridge, especially at the two vertical edges of the 
ZGNR bridge (see Fig. 0(b) and (c)). The spatial distri- 
butions of the quasi-bound states show the antibonding 
(Fig. [6] (b)) and bonding (Fig. [6](c)) states correspond- 
ing to the resonant tunneling peaks B and C marked in 
Fig. EH(a). This bonding and antibonding feature also ex- 
plain why the positions of the resonant tunneling peaks 
B (C) decrease (increase) as the length Ntc increases. 

Finally, we investigate the effect of the widths of the 
right lead on the transport property. We fix the width 
of the left lead at Nl = 20 while tuning the width of the 
right lead. Figure [7l (a) shows the contour plot of the 
conductance as a function of the Fermi energy Ef and 
the length of the ZGNR bridge N T c for the width of the 
right lead Nr = 10. The figure shows that the first kind 
of tunneling peaks disappears but the second kind of tun- 
neling peaks still can be observed, because the latter is a 
consequence of tunneling through the quasi-bound states 
localized in the middle bridge. The effect of the width of 
the lead on the second kind of tunneling processes is much 
weaker than on the first kind of tunneling processes. For 
the wider right lead Nr = 40 in Fig. [71 (b), the situation 
is similar with Nr = 10. But when the Fermi energy is 
larger than the bottom of the second subband of the right 
lead, a strong tunneling process appears at Ef ~ 0.08£. 
With a wider or narrower right lead, the second kind of 
tunneling processes display similar behavior as the length 
(Ntc) of the bridge increases. The tunneling processes 
induced by the quasi-bound antibonding states results 
in tunneling peaks with an energy that is approximately 
0.04t higher. Those peaks shift to lower energy as the 
length of the bridge Ntc increases. In contrast, opposite 
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FIG. 7: (Color online) The same as Fig [5] (a), but (a) and (b) 
are for Nr=10 and 20, respectively. 

behavior as compared to those induced by the antibond- 
ing states, i.e. the blue shift. 



tion of the Fermi energy the tunneling current can be 
turned on and off periodically making this device suitable 
for transistor action. We demonstrated that resonant 
tunneling can be realized utilizing the geometrical effect 
of such bended GNR structures even in the absence of 
any nanostructured gates. The strong resonant tunneling 
processes arise from the quasi-bound states that localize 
at the middle part and can be clearly seen in the con- 
ductance spectra. The spectra exhibits rich structures 
coming from the different kind of quasi-bound states in 
the middle part. The resonant tunneling processes can be 
tuned by changing the Fermi energy and the width of the 
left and right leads. Our theoretical results demonstrate 
that resonant tunneling can be realized utilizing the ge- 
ometry of the graphene nanoribbon without any nanos- 
tructured electric gates. We believe that our theoretical 
results are interesting both for the basic physics and for 
the potential application of electronic devices based on 
graphene nanoribbons. 



IV. CONCLUSIONS 
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